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Abstract. In grand unified theories (GUT), non-universal boundary conditions for the gaugino masses 
may arise at the unification scale, and affect the observability of the neutral MSSM Higgs bosons (h/H/A) 
at the LHC. The implications of such non-universal gaugino masses are investigated for the Higgs boson 
production in the SUSY cascade decay chain g — > qq, q — > x~2q, ill — > x~ih/H/A, h/H/A —* bb produced in 
pp interactions. In the singlet representation with universal gaugino masses only the light Higgs boson can 
be produced in this cascade with the parameter region of interest for us, while with non-universal gaugino 
masses heavy neutral MSSM Higgs boson production may dominate. The allowed parameter space in 
the light of the WMAP constraints on the cold dark matter relic density is investigated in the above 
scenarios for gaugino mass parameters. We also demonstrate that combination of representations can give 
the required amount of dark matter in any point of the parameter space. In the non-universal case we show 
that heavy Higgs bosons can be detected in the studied cascade in parameter regions with the WMAP 
preferred neutralino relic density. 

PACS. 12.60.Jv Supersymmetric models - 95.35.+d Dark matter - 14.80.Cp Non-standard-model Higgs 
bosons - 14.80.Ly Supersymmetric partners of known particles - 12.10.Kt Unification of couplings; mass 
relations 



>■ ! 1 Introduction 



Most studies concerning the discovery potential for find- 
ing supersymmetric Higgs bosons at the LHC, have been 
done in models with universal gaugino masses at the unifi- 
cation scale. There is, however, no compelling theoretical 
reason for such a choice. In grand unified supersymmet- 
ric models, which include an SU(5) grand unified model, 
non-universal gaugino masses are generated by a nonsin- 
glet chiral superfield <P n that appears linearly in the gauge 
kinetic function /(<£). The function f{<F) is an analytic 
function of the chiral superfields in the theory [1]. It 
should be noted that the chiral superfields <P consist of 
a set of gauge singlet superfields <P and gauge nonsin- 
glet superfields <P n , respectively, under the grand unified 
group. If the auxiliary part F$ of a chiral superfield <P in 
the f(<P) gets a VEV, then gaugino masses arise from the 



coupling of with the field strength superfield W a . 

The Lagrangian for the coupling of gauge kinetic function 
with the gauge field strength is written as 



C 9 k = / d 2 6f ab ($)W a W b + H.c. 



(1) 



where a and b are gauge group indices [for example, a,b = 
1,2, ...,24 for SU(5)], and repeated indices are summed 
over. The gauge kinetic function / Q &(<£) is 



(f,n 



(2) 



where as described above the <P S and <P n are the singlet 
and nonsinglet chiral superfields, respectively. Here /o(^ s ) 
and / n (^ s ) are functions of gauge singlet superfields <P S , 
and Mp is some large scale. When F<p gets a VEV (F&), 
the interaction (p} gives rise to gaugino masses: 
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>-\ a \ b + H.c. 



(3) 



where A a,fc are gaugino fields. It should be noted that the 
U(l), SU(2), and SU(3) gauginos are denoted by Ai, A2, 
and A3, respectively. 
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Since the gauginos belong to the adjoint representa- 
tion of the gauge group, in the case of SU(5) for example, 
<P and F$ can belong to any of the following representa- 
tions appearing in the symmetric product of the two 24 
dimensional representations of SU(5): 

(24®24) Sj , mm = 10 24 75 9 200. (4) 

In the minimal case (which is the simplest one too), ^ 
and F$ are assumed to be in the singlet representation of 
SU(5). This corresponds to equal gaugino masses at the 
GUT scale. However, <P can belong to any of the nonsinglet 
representations 24, 75, and 200 of SU(5). In that case, the 
gaugino masses are unequal but related to one another via 
the representation invariants [2113111]. It should be kept in 
mind that an arbitrary combination of these different rep- 
resentations is also allowed. We shall first study the case of 
each representation separately, and then the effect of com- 
bining two of these. In Table Q] we display the ratios of re- 
sulting gaugino masses at tree level as they arise when 
belongs to various representations of SU(5). Clearly, the 
nonsinglet representations have characteristic mass rela- 
tionships for the gauginos at the GUT scale. The resulting 
relations at the electroweak scale, using the renormaliza- 
tion group (RG) evolution at the one-loop level are also 
displayed. 



Table 1: Ratios of gaugino mass parameters at the GUT 
scale in the normalization M3(GUT) = 1 and at the elec- 
troweak scale in the normalization Ms(EW) = 1 at the 
one-loop level. 
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The phenomenology of supersymmetric models with 
nonuniversal gaugino masses has been considered e.g. in 
P61I71I8PIT01ITT1IT21. Studies were performed also in the 
context of supersymmetric dark matter [13. 14, 15, 16p7] 
fT5UT5] . The phenomenology of supersymmetric models de- 
pends crucially on the compositions of neutralinos and 
charginos. Hence, it is extremely important to investi- 
gate the changes in the experimental signatures with the 
changes in the composition of neutralinos and charginos 
which arise because of the non-universal gaugino masses 
at the GUT scale. In this study our focus will be on the 
experimental signatures of neutral Higgs bosons appear- 
ing in the cascade decay chains of the squark and gluinos 
produced at the LHC. If the squarks and gluinos are light 
enough, their production cross sections are large at the 
LHC. The light neutralinos xi,2 are typical decay prod- 
ucts of g and q. The neutral Higgs bosons can be produced 
in the decay of xi , if the mass difference between xi and 
Xi is large enough. As the production rate is largely inde- 
pendent of the value of tan/3, these production channels 



have been found particularly interesting at the LHC to 
cover the difficult region of low and medium tan/3 values 
20, 21j . Production of Higgs via xi — > Xih/H/A in models 
where gauginos are in singlet and 24 representation were 
studied in [5]. Also, recently the role of gaugino mass non- 
universality (in a model independent framework) in the 
context of Higgs productions has been discussed in [22\. 
It was found for sample parameters that only the light- 
est Higgs could be produced in the model with singlet, 
while also the heavy Higgses, H and A, could be produced 
from the cascade in the model with 24 representation. 
The gaugino mass parameters indeed indicate that it is 
more probable to have the cascade decay including heavy 
Higgses in the model with gauginos in 24. Namely, if the 
fi parameter, mixing Higgs superfields in the MSSM su- 
perpotential, is large, the lightest neutralino using either 
representations is close to bino, while the second lightest 
is close to wino. From Table 1, in the model with singlet 
the mass difference of the neutralinos is then close to the 
mass of the lightest neutralino, since m X2 ~ 2m Xl , while 
in the model with 24, m X2 ~ 6m Xl and the mass differ- 
ence thus five times the mass of the lightest neutralino. 
Note that the cascade is not possible with small values of 
the //-parameter, since in such a case %2 and xi are similar 
in mass. From Table 1, one can also note that the cascade 
in the cases of representation 75 is not probable. 

In this work, the Higgs boson decay to b quarks, 
H/A/h — > bb, with a large branching fraction (~0.9) is 
selected. In the multi-jet events this decay mode suffers 
from a large combinatorial background, but structuring 
the event into two hemispheres according to the direc- 
tions of the decaying gluinos and squarks from the hard 
process, this background can be efficiently reduced. Pro- 
duction of the light Higgs boson h in SUSY cascades with 
bb final state has been studied in 23 in the CMS de- 
tector with the assumption of universal gaugino masses. 
The event selection is based on a requirement of four en- 
ergetic jets, large missing Et, separation of the jets into 
two hemispheres and the reconstruction of the Higgs bo- 
son mass from two jets tagged as b jets. Similar analysis 
method is used here. The main phenomenological differ- 
ences between approaches of Ref. |23j and this work are 
the different branching fractions of the SUSY particles. 
In the singlet representation only the light Higgs boson 
h can be produced whereas in the representation 24 the 
X2 — ► H/Axi decays dominate in most of the parameter 
space. Therefore considerable differences can be expected 
in the final mass distribution as well as the obtainable 
significance of the signal. In Ref. [5] the 4b final state 
from H/A/h — > bb decay from both of the cascades has 
been proposed as the search channel. This mode however 
is subject to a large combinatorial background from bb 
combinations due to real or fake b jets in the cascade not 
containing the Higgs boson. 

Combining information from different sources is essen- 
tial for utilizing future facilities in best possible way. In 
addition to the laboratory studies, relevant input is ob- 
tained from the dark matter searches: the WMAP satellite 
has put precise limits on the relic density. Supersymmetric 
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theories which preserve R parity contain a natural candi- 
date for the cold dark matter particle. If the lightest neu- 
tralino is the lightest supersymmetric particle (LSP), it 
can provide the appropriate relic density. In many SUGRA 
type models the lightest neutralino is bino-like, which of- 
ten leads to too large thermal relic density, as compared to 
the limits provided by the WMAP experiment [24, 25J. Of 
the possible representations of SU(5) GUT all but the sin- 
glet representation produce nonuniversal gaugino masses. 
Thus the resulting neutralino composition changes from 
the usual universal gaugino mass case [5]- Depending on 
the representation where the gauginos are, this leads to 
different regions in the parameter space where WMAP 
constraints are satisfied. It is also shown that by combin- 
ing representations, observationally allowed dark matter 
parameter regions are not very restricted, if parameters 
are suitably chosen. 

Here we are interested in the case when both the relic 
density constraints from WMAP can be satisfied and the 
heavy Higgses can be produced in the cascade. An im- 
portant contribution to the neutralino annihilation comes 
from the enhancement due to the Higgs resonance. If we 
denote the mass of the heavy Higgs, either H or A, by tua, 
the enhanced annihilation occurs when = 0.5mA- 

Since the mass difference between the two lightest neu- 
tralinos is 0.5mA in the model with singlet and 2.5mA in 
the model with 24, it is obvious that this mechanism is 
available in the latter case, but not in the case with sin- 
glet, if we want to study the Higgs production via cascade 
at the same time. Indeed, we find easily parameter sets 
with both neutralino dark matter and the cascade to all 
the Higgses h, H, and A in the model with 24 but not in 
the model with singlet. Thus in this work we are especially 
interested in the model with representation 24, for which 
also the production cross sections for heavy Higgses were 
found reasonably large in [5]. In order to have both rele- 
vant region for dark matter and to produce heavy Higgses 
in the cascade decays, the model with 24 is special among 
those with nonuniversal SU(5) representations. We shall 
see later that in the model with 200, the relic density is 
never large enough to produce the dark matter solely. In 
75 the mass difference between the two lightest neutrali- 
nos is too small for the cascade decay in the parameter 
region with relic density in the WMAP region. 

In the following we first discuss in section 2 the pa- 
rameter regions for observationally acceptable dark mat- 
ter in different SU(5) representations. In section 3 we con- 
centrate on representation 24 and specifically study the 
parameter region, where the Higgs production from neu- 
tralino decay is most promising, and the signal and back- 
ground processes are discussed. The detector simulation 
is shortly presented in Section 4. The event selection and 
the signal to background ratio are discusssed in Section 5, 
and the background determination is elaborated in Section 
6. Section 7 describes the systematic uncertainties on the 
background determination. Discovery potential is given in 
Section 8 and the conclusions in Section 9. 



2 Dark matter in SU(5) representations 

We start the description of the relic density for four dif- 
ferent SU(5) representations with the representation 24, 
which is most interesting for us here, as explained in the 
introduction. We then shortly discuss the representations 
75 and 200, and compare with the singlet, which has been 
studied extensively in the literature. Then we discuss com- 
bination of singlet representation with 200. 

The amount of thermal relic density in the represen- 
tation 24 is presented in Fig. Q] for four different sets 
of GUT scale parameters. The spectrum was calculated 
with SOFTSUSY 26J and relic densities and constraints with 
micrOMEGAs (2"7 ] l28 p 29j . The dark shaded areas represent 
larger relic density than the lighter areas. For the relic 
density, we use here the WMAP combined three year lim- 
its HUES] 



^CDAlh 2 



iin c U+ 00976 



(2a). 



(5) 



For the particle masses, the following limits are applied 
[29] : me R > 99.4 or 100.5 GeV depending if the lightest 
neutralino mass is below or above 40 GeV, m^ R > 95 
GeV, TOf x > 80.5 to 88 GeV depending on the lightest 
neutralino mass (from 10 to 75 GeV), mp i > 43 GeV, 
and m^± > 73.1 to 103 GeV depending on the sneutrino 
masses (from 45 to 425 GeV). The curve m^ = 114 GeV 
is depicted in the figure. For the shown parameter region, 
when otherwise experimentally allowed, Higgs is always 
heavier than 91 GeV, which is the Higgs mass limit in 
MSSM for tan/3 > 10 assuming maximal top mixing [30J. 

The most striking feature in Fig. [T] is the valley of the 
low relic density area around Mi ~ —300 GeV. There 
the annihilation of neutralinos proceeds to quarks and the 
minimum occurs at the Z peak providing an efficient anni- 
hilation. Outside of the valley the relic density rises, over- 
closing the universe. In the representation 24 the lightest 
neutralino is very bino-like, and the WMAP preferred re- 
gion tends to be quite narrow. 

If the overall relic density is low, the crucial differ- 
ences that control the preferred regions of the parameter 
space come mostly from the composition of the lightest 
neutralino. An important single aspect that separates the 
different parameter sets in Fig. [1] is the presence of the 
higgsino component in the lightest neutralino, which in- 
creases annihilation into the gauge boson pairs. Increasing 
the value of tan (3 decreases the value of the supersymmet- 
ric higgsino mass parameter fj,, as can be seen from the 
electroweak symmetry breaking (EWSB) condition 



H 2 = 



'H, 



tan 2 (3 1 



tan^ (3-1 



' z- 



(G) 



where tan/3 = (H^)/(H^) = v u /vd- This results in a big- 
ger higgsino component in the lightest neutralino, which 
then annihilates more efficiently. The effect is seen in com- 
paring Figures [Taland fTbl Since the relic density in Fig.llbl 
is smaller than in Fig. [Ta| in general, also annihilation to 
the lightest Higgs suffices to lower the relic density to the 
WMAP preferred region. Thus the low relic density re- 
gion is wider than in Fig. [Taj with bands corresponding 
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Fig. 1: Relic density f2 x h 2 in the representation 24 in (M 2 ,mo) plane. The dark shaded areas represent the larger 
relic density, wmap denoted filling is the WMAP preferred region, lep shows an area next to the axes, where the 
experimental mass limits are not met, rge shows an area where there is no radiative EWSB, and Isp the area where 
neutralino is not the LSP. h gives the line of m/j = 114 GcV. The area between two bsg contours is allowed by b — » s-f 
limits in (a), (b), and (d), and above the contour in (c). 



to both Z and h. The higher tan (3 also increases the area 
at the low too, where the radiative EWSB is not possible 
(marked rge, below the area where neutralino is not the 
LSP).Q 

The same mutual sign for the [i parameter and the 
gaugino mass M\ gives a slightly lower mass for the light- 
est neutralino than in the case with opposite sign. Since 
in the representation 24 the signs of M\ and M 2 are 
the same, this effect can be seen comparing figures Hal 
and [Tel Thermal relic density is proportional to the mass 
of the particle. The b — > sj constraint is also promi- 
nent in the negative \i case [31,32,33,34,35:. We have 

1 It can be noted that the Higgs funnel region, which is 
present in singlet case, does not exist in 24 representation at 
large tan (3 and large mo values. This is due to the neutralinos, 
which are too light to annihilate to the heavy Higgs bosons at 
the resonance. 



used the two sigma world average of BR(b — > sj) = 
(355 ± 24+% ± 3) x 1(T 6 for the branching fraction [36]. 
For example in Fig. [Tc] the visible area is mostly excluded 
by the 2a limit. In other figures the b — > sj allowed region 
is between the frsg-denoted contours. 

Switching on the trilinear scalar coupling parameter 
Aq increases the value of \i through the RG evolution of 
the squared soft Higgs mass parameters, therefore reduc- 
ing the higgsino component in the lightest neutralino. This 
can be seen comparing Figures [Tc] and [id] at the low | M 2 \ 
values. In the Fig. [Tc] the higgsino component of the LSP 
is larger at the low jil^l-high mo region, so the resulting 
relic density is lower than in Fig. lid] where the lightest 
neutralino is almost purely binoo 

2 In Fig. Ildl the relic density can rise up to the values of 80; 
the shading reaches maximum saturation already in 50 in order 
to ease comparison with other figures. 
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Fig. 2: Relic density Q x h 2 in the representation 75 in 
(M2, too) plane. Explanation for the key is as in Fig.[TJ The 
parameters used here are tan/3 = 10, sgn(/x) = — 1, Ao = 
1 TeV. 



In Fig. [2] the contours of preferred thermal relic den- 
sity in the representation 75 are plotted for one set of 
(GUT scale) parameters, for which appreciable region for 
dark matter within the WMAP limits is found, and which 
are not forbidden by the b — > 57 constraints. Since the 
higgsino component in the representation 75 is large [5], 
the resulting relic density is low, and most of the param- 
eter space is not overdosed by the WMAP limits. Also 
the co-annihilations with the lightest chargino reduces the 
relic density, since the lightest neutralino and chargino are 
nearly mass degenerate in the higher Mi part of the pa- 
rameter space. This is also seen in the Fig. [2] at high too, 
where the lightest chargino becomes the LSP for specific 
Mi values. In the low M 2 region the EWSB condition 
pushes the \i high, which in turn decreases the higgsino 
component in the lightest neutralino making it mostly a 
bino. The lightest neutralino and chargino are not degen- 
erate anymore, and the relic density increases in low Mi 
area, which enables the emergence of the WMAP preferred 
region in the parameter space. The second lightest neu- 
tralino can annihilate also directly into gauge bosons in 
this parameter region. At the low mo region the LSP can 
be the lighter stop. Again, increasing tan /3 enhances the 
higgsino component leading to lower relic densities in gen- 
eral. 



In Fig.[3]the contours of preferred thermal relic density 
in the representation 1 are plotted for four different sets of 
(GUT scale) parameters for the reference. For the chosen 
parameters the WMAP preferred regions are found near 
the Mi (i.e. m 1/2 in mSUGRA and CMSSM language) 
and too axes0 The b — > sj limits exclude large area in 
the figures [3c] and I3d( where the sign of the ^-parameter 
is negative. 

In the representation 200 the higgsino mixing is large 
5J and importantly, the bino-component is very small. 
Also the co-annihilations with the lightest chargino, which 
is nearly mass degenerate with the lightest neutralino, re- 
duce the relic density substantially, and the resulting relic 
density is tiny. In contrast to the 75 dimensional case, the 
/i parameter decreases with decreasing Mi, so the bino 
component does not get very large. The b — > sj con- 
straint is largely within acceptable limits (for negative \x 
a bit larger M 2 and too are required than for positive /i) 
and the LSP is mostly the lightest neutralino. Typically 
the Higgs mass is well above 91 GeV limit for tan/3 >, 10, 
and the 114 GeV line is around M 2 ~ 600 GeV. For the 
lower values of tan /3, the Mi parameter must be above 
1 TeV. The parameter space suitable for finding partial 
neutralino dark matter can be extended both in too and 
in Mi beyond 1 TeV for both signs of n, but neutralino in 
this representation can never be the only source of dark 
matter. 

It is possible also to have many representations simul- 
taneously. For example, the Fig. [4] shows the relic density 
when two representations, 1 and 200, are combined. If one 
allows the gaugino masses to be formed from two repre- 
sentations simultaneously, then two of the gaugino masses 
can be chosen freely, and only one is fixed by the repre- 
sentations. In Fig. [4] the parameter M3 is chosen to have 
a specific value (M 3 = 700 GeV in Fig. \M and M 3 = 400 
GeV in Fig.l4bp and Mi is allowed to vary. The overall relic 
density is very low, as implied by the representation 200. 
However, there is a narrow peak where the relic density 
is enormous. This is the area, where the gaugino masses 
combine in such a way that the lightest neutralino is ex- 
tremely light. Therefore the annihilations are very much 
suppressed leading to high relic density. The position of 
the area where this happens depends on the representa- 
tions: e.g. for the representation 1 200 the parameter 
Mi can be expressed as 

Mi = m.(i-iQHW0-»> = , m ,-m,„ (7) 

where the numbers come from the representation invari- 
ants. Therefore the minimum of the lightest neutralino 
mass is positioned around M2/M3 ~ 8/9. The b — > 57 
constraint is within acceptable limits and the Higgs mass 
is above limits (in Fig. [4b] m/,, > 113.5 GeV in otherwise 

3 The hyperbolic branch/focus point region [371138]. which 
requires a substantially larger mo, is not found for this setup 
(moreover, the low scale parameterization might be more rel- 
evant in describing this region anyway [39]). Also, the A- 
annihilation funnel is absent, since it requires larger tan/? ~ 
45 - 60. 
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Fig. 3: Relic density f2 x h? in the representation 1 in (M 2 ,mo) plane. Explanation for the key is as in Fig.[TJ 



allowed regions) . By increasing the value of the M3 mass 
parameter one effectively increases the stop masses (due 
to the RG evolution), which then increases the lightest 
Higgs mass. For example, the rrih > 117 GeV always in 
the Fig. Ha) Chargino being the LSP excludes large areas 
in Fig.[3al The smaller area around M 2 ~ 800 GeV at the 
low mo region is due to the stau LSP. The A-threshold is 
also shown as a small ring of higher relic density around 
M 2 ~ 770 GeV. 



Because with two representations there are two free 
gaugino masses, the light LSP region is always present. 
Therefore, if the relic density is low in general, there must 
also be a stripe of WMAP preferred region available. Chang- 
ing the M3 mass value the WMAP preferred stripe can 
be moved through the parameter space, thus providing a 
method to fulfill the WMAP constraints. 



3 Higgs in the neutralino cascade and dark 
matter 

Having discussed the basic features of the SU(5) non- 
universal gaugino mass scenario, and the region of param- 
eter space preferred by the thermal relic density, we de- 
scribe the numerical features of the cross section of the sig- 
nal events in the representation 24. In Fig.[5]the cross sec- 
tion for g -> qq, g -> x~2<7, X~2 -> Xih/H/A, h/H/A -> bb 
produced in the pp interactions in the (to a, tan j3) (?ig- 
ures EU [5cl [5e]) and (m^ra^) (figures l5b" l l5dl l5fjl planes 
are plotted. The solid (green) fill denotes the WMAP 
preferred relic density region. Also the Higgs 114 GeV 
mass contour is plotted with the contours of constant 
cross section (larger tua values correspond to larger mh)- 
The values of the parameters are (fi — +700 GeV/c 2 , 
TOq = 600 GeV/c 2 , m i = 350 GeV/c 2 , and the trilin- 
ear coupling for the top sector is chosen to be A t = 800 
GeV). In the (m^, tan/3) figures the gluino mass is cho- 
sen as TOg = 770 GeV/c 2 , and in the (to^jTOq) plane 
tan/3 = 10. In the (TO^jtan/3) plane the cross section 
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is largest for the H Higgs production. Also the stripe of 
WMAP preferred relic density passes through the large 
cross section area. For the A production, the cross sec- 
tion is somewhat smaller, while for the h production it is 
substantially smaller for the WMAP preferred area. The 
near- horizontal WMAP stripe around tan/3 = 8 (begin- 
ning at tiia — 130 GeV) corresponds to the closing of the 
Higgs resonance in the LSP annihilation: in the region be- 
low that line the lightest neutralino mass is more than half 
of the light Higgs boson mass, and thus the annihilation 
never occurs at the resonance. Heavier A correspond to the 
larger relic density. The relic density is lowest in between 
the two relic density stripes, the minimum occurring at 
the A-peak. 

In the (m J 4,m§) plane the large H production cross 
section is again passed by the preferred relic density stripe. 
The Higgs 114 GeV limit divides the parameter space in 
two (vertical line in the figures EH [HH1 H3 it should be 
noted, though, that for the tan/3 = 10 used here the ac- 
tual Higgs boson mass limit can be as low as around 90 
GeV). Also here the relic density is the lowest in between 
the two relic density stripes. The horizontal kink in the 
relic density stripe around rrig ~ 770 GeV corresponds 
to the closing of the Higgs resonance in the LSP anni- 
hilation: in contrast to the (m^, tan/3)-plots, the Higgs 
resonance is open below that line due to the decreasing 
of the lightest neutralino mass. The A and h cross sec- 
tions follow the pattern of the (m a, taxi j3) plane, the A 
production cross section being slightly smaller than for H 
while h production is the one with the lowest cross section. 
The cross sections and relic densities were calculated with 
PR0SPIN0 @0], SDECAY [41], HDECAY [42] and SuSpect [43] 
and mi cr OMEGAs [2"M2"5ll2T] . 

Table [2] shows the mass values for \i an d Xi an d the 
branching fractions for the xi ~> h/H/Axi decay channels 
in the representations 1 and 24 with hia = 190 GeV/c 2 
and in the representation 24 with niA = 210 GeV/c 2 



with tan/3 = 10. For the 24-dimcnsional representation 
all Higgs boson channels are possible. The branching frac- 
tion is largest for the CP-even neutral Higgs boson H. 
As the gluino mass increases, the branching fraction to 
X2 — > h/H/A — ► xi decreases. The signal cross sections are 
obtained applying the % 2 — > h/H/Axi and h/H/A — > bb 
branching fractions on the total SUSY production cross 
section in the chosen parameter point. In the singlet case, 
the Trig = 770 GeV is too small to allow kinematically the 
cascade decay. We take rrig = 850 GeV, with which the 
decay through the light Higgs boson h is possible. 

4 Event generation 
4.1 SUSY cascades 

The event generation and simulation were performed and 
the results were made public in the CMS framework |44j . 
The Higgs boson production in SUSY cascades was gener- 
ated with PYTHIA [45] using a squark and gluino produc- 
tion with general MSSM simulation. The CTEQ5L struc- 
ture functions |46] were used. The simulated point was 
H = +700 GeV/c 2 , m h = 600 GeV/c 2 , = 350 GeV/c 2 , 
and the trilinear coupling for the top sector was chosen 
to be A t — 800 GeV. The gluino mass was chosen as 
m g = 770 GeV/c 2 . The U(l), SU(2) and SU(3) gaug- 
ino masses were chosen separately for each scenario, val- 
ues shown in Table [3] PYTHIA version 6.325 was used 
within the CMKIN framework [47]. No forcing of decay 
channels was used. After the event generation, the events 
were sorted by the Higgs boson content, giving a sample 
of signal events, and a sample of SUSY background events 
with no Higgs bosons in the event. The SUSY production 
of multi-jet events is expected to be the main source of 
background. 

In the Singlet representation a cross section of 18.3 
pb was used for the SUSY cascade (including Higgs boson 
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Fig. 5: Contours of constant cross section for pp — > h/H/A 
(niA,mg) (right) planes. The solid fill denotes the WMAP p 
shows the Higgs 114 GeV mass limit (higher vha values corr< 
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referred relic density region and the dense dashed curve 
jspond to a higher m^). 
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Table 2: Mass values for xi and X2 and the branching fractions for the \2 ~ * h/H/A/fi decay channels in the 
representations 1 (nig — 850 GeV) and 24 (trig = 770 GeV) with hia = 190 GeV/c 2 and in the representation 24 
with m A = 210 GeV/c 2 and tan/3 = 10. 



Representation, hia 


m* (GeV/c^) 


m x - 2 (GeV/c^) 


BR(x2 hxi) 


BR(x2 -> Hxi) 


BR(x2 -» Axi) 


1, 190 GeV/c 2 


120.6 


246.8 


0.897 








24, 190 GeV/c 2 
24, 210 GeV/c 2 


50.5 


307.4 


0.00478 


0.718 


0.205 


50.5 


307.5 


0.00144 


0.739 


0.161 



production) with niA = 190 GeV/c 2 , tan/3 = 10. The Higgs 
boson production cross section is 4.13 pb for this point. In 
the 24-dimensional representation with niA = 190 GeV/c 2 , 
tan/3 = 10, the cross sections are 23.1 pb and 5.07 pb, and 
with m A = 210 GeV/c 2 , tan/3 = 10, 23.1 pb and 4.97 pb, 
respectively. 



4.2 Standard Model backgrounds 

The tt, QCD multi-jet, Z+jets and W+jets production are 
the main SM backgrounds for the Higgs boson searches 
with H/A/h — » bb in SUSY cascades. The tt events were 
generated with PYTHIA [45]. The production cross sec- 
tion was normalized to the NLO cross section of 840 pb. 
The QCD multi-jet background for the light Higgs boson 
h in the SUSY cascades has been studied with full sim- 
ulation in Ref. [23] and has been found to be negligible 
requiring energetic jets, two b jets and large missing Et 
(> 200 GeV). This background particularly sensitive to 
the detector performance through the jet measurements. 
Therefore it was studied in this work with the new full 
simulation [35] for the E™ lss cut. Events for the full sim- 
ulation study were generated with PYTHIA in two bins: 
120 < p T < 170 GeV/c and 170 < p T < 230 GeV/c. The 
Z+jets and W+jets events can contribute to the back- 
ground through large missing Et from Z — > vv decays 
and from lost leptons. These backgrounds have been in- 
vestigated for the inclusive SUSY searches with ALPGEN 
[39] and PYTHIA simulation. The Z+jets and W+jets 
backgrounds have been found to be roughly on the same 
level as the tt background after the basic selections for 
the inclusive search: E™ ISS > 200 GeV and four jets with 
Et > 30 GeV [50]. The requirement of two b jets with 
Et > 100 GeV will suppress these backgrounds well be- 
low the tt background and were therefore not simulated 
in this work. 



5 Detector simulation 

The response of the CMS detector was simulated with 
the fast simulation package FAMOS [51] . The pile-up cor- 
responding to the luminosity of 2x 10 33 cm _2 s~ 1 was in- 
cluded. The physics objects were reconstructed with the 
standard methods available in the CMS reconstruction 
software. The primary vertex was reconstructed and se- 
lected with an algorithm searching for the highest sum of 
the transverse momenta of the associated tracks [52! . Jets 



were reconstructed in a cone of 0.5. The b tagging was 
performed with a robust track counting based B- Tagging 
Algorithm of CMS [53]. The QCD MET reconstruction 
and the MET cut efficiency was studied in addition with 
CMS full simulation [48]. 



6 Event selection 

The signal events are characterized by many energetic 
hadronic jets. Due to the mg > rriq condition, the gluino 
decays to a quark and the corresponding squark. There- 
fore the cascade producing the Higgs boson contains two 
energetic jets from the decays of the gluino and the squark 
and two b jets from the decay of the Higgs boson. At least 
two jets are produced in the opposite cascade. The two 
LSP's (xi) at the end of the two cascades lead to large 
missing Et- These characteristics can be used to suppress 
the SM backgrounds, for which the jet multiplicities and 
missing Et are typically significantly smaller. 

6.1 Trigger 

The events are supposed to be triggered with a multi-jet 
trigger at the Level-1 (LI) and with a combined multi-jet 
plus missing Et trigger at the High Level Trigger (HLT) 
54.55]. The Et thresholds for the four jets at the trigger 
level were taken to be 80, 50, 30 and 10 GeV. The thresh- 
old for the missing Et calculated from the calorimeter 
cells was fixed to 135 GeV. The Level-1 trigger was not 
simulated. The effect of this choice is negligible since the 
efficiency of the multi-jet trigger on the SUSY events has 
been found to be close to 100% 23J. In this work the HLT 
efficiency in the SUSY events is found to be 67%. The 
tt background is efficiently reduced already at the trigger 
level with the missing Et cut. The HLT efficiency for this 
sample was found to be 3.8%. 



6.2 Jet selection 

Figure [6] shows the Et distributions for the three most 
energetic jets in the signal, in the SUSY background and 
in the tt events. For the off-line analysis, events with four 
jets with Et > 50 GeV were first searched for. At this 
level, the E T cuts were set to 250, 200, 100 GeV for the 
three most energetic jets. These cuts are more stringent 
than those used in Ref. [23] and thus can help to suppress 
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Table 3: Gaugino mass parameters for the studied points. 



Representation, niA 


U(l)(GeV/c 2 ) 


SU(2)(GeV/c 2 ) 


SU(3)(GeV/c 2 ) 


1, 190 GeV/c 2 
24, 190 GeV/c 2 
24, 210 GeV/c 2 


122.9 
-44.1 
-45.54 


245.3 
-264 
-272.7 


850 
610 
630 



more efficiently the QCD multi-jet and W+jets and Z+jets 
backgrounds. The tt background is reduced to 35% with 
the E^* cuts. Two b jets are expected from the Higgs 
boson decay. Figure [7] shows the px distributions for the 
two b quarks from H/A — > bb. Efficiency for tagging one 
selected jet as a b jet in the signal events, with at least 
one Higgs boson decaying to b quarks, was found to be 
about 30%. In addition to the possibility of mistagging 
a hadronic jet in the SUSY background, there are several 
sources of genuine b jets, like %2 — > Z^i, Z — > bb and gluon 
splitting processes in hadronic jets. As a consequence, the 
tagging efficiency in the SUSY background was found to 
be substantial, about 13%. 

Figure [5] shows the distribution of the missing Et for 
the signal, for the SUSY background and for the tt events. 
Large missing Et is typical for SUSY events due to the 
presence of neutralinos in the final state. The tt back- 
ground is further suppressed with this cut by a factor of 
~ 5. The QCD multi-jet production contributes to the 
background through missing Et due to detector resolu- 
tion and mis-measurement of jets. The missing Et se- 
lection was studied with the new full simulation pack- 
age gS]. Efficiency for E™ 88 > 200 GeV was found to 
be 4xl0~ 5 for 120 < p T < 170 GeV/c and 9xl0~ 5 for 
170 < pt < 230 GeV/c. No trigger simulation was in- 
cluded in the estimation of this background. 

Due to the large jet multiplicity and the mistagging 
possibility, the combinatorial background from false jet 
assignment can be expected to be large for the recon- 
struction of the Higgs boson invariant mass from the b 
jets. This combinatorial background can be significantly 
reduced with a method called a hemisphere separation 
technique, devised in |56j . This method separates the jets 
into two "hemispheres" where the initial gluino or squark 
jet originates, starting from the most energetic jet in the 
two hemispheres. The assignment of a jet to the axis is 
performed exploiting the angles and invariant masses. The 
method has been shown to be effective in removing the 
combinatorial background arising from jet multiplicity in 
SUSY events in Ref. [13]. 

After the separation to hemispheres the b jets were 
selected in each hemisphere separately. Due to the large 
mass of the decaying neutralino, the Higgs boson, in par- 
ticular the light Higgs boson h, receives a large boost lead- 
ing to a small separation angle between the b jets. The 
AR (AR = ^/(Ary) 2 + (A^) 2 ) between the two b quarks 
from the Higgs boson, shown in Fig. [5] for the light Higgs 
boson, is peaked around Z\R ~1. Therefore the pair of b 
jets with the lowest Z\R value was selected in each hemi- 
sphere. Among these two combinations the one with lower 
AR value was taken to present the Higgs boson. 



Figure [10] shows the bb mass distribution for the sig- 
nal in the singlet representation with hia = 190 GeV/c 2 . 
Only the light Higgs boson can be produced due to the 
small mass difference between the two neutralinos. The 
mass distribution with correct b-jet mass assignments is 
also shown in the figure. The tagged b jet is defined as a 
genuine b jet if a b quark is found within AR < 0.4 from 
the jet axis. The fraction of correct jet assignments is 29%. 
Large fraction of the false b-jet assignments is distributed 
to large mass values. A tail at large masses is visible also 
in the mass distribution with correct assignments. This 
tail is due to overlapping jets in the multi-jet environ- 
ment. Figure [TT] shows the bb mass distribution for the 
signal in the representation 24 with hia = 190 GeV/c 2 . 
For this parameter point the mass difference between the 
neutralinos is large and the decay branching fractions to 
heavy Higgs bosons are significant. The fraction of the 
light Higgs boson after all selection cuts is 3.5%. For the 
heavy Higgs bosons, better signal mass resolution and sig- 
nal visibility was obtained with a significantly harder cut 
for the selected b jets, Et > 100 GeV, as is shown in 
the Fig. [TT] About 87% of the signal is due to the heavy 
CP-even Higgs boson H. The mass distribution with cor- 
rect b-jet mass assignments and with Et > 100 GeV is 
also shown in the figure. The opening angle in the Higgs 
boson decay is larger than for the light Higgs boson h 
in the singlet representation making more favorable the 
jet reconstruction and b-jet identification. The fraction of 
correct jet assignments was found to be 39%. 

6.3 Signal over background ratios 

Figure [TU shows the bb mass distribution for the signal 
in the singlet representation with hia = 190 GeV/c 2 and 
tan/3 = 10 and for the total background. The tt back- 
ground is shown separately in the figure. Setting an upper 
bound of AR < 1.5 could reduce this tail as can be seen 
from Fig. [13] To further reduce the backgrounds corre- 
lations between the reconstructed bb pair, presenting the 
Higgs boson, and the missing Et was studied. The missing 
Et originates from the \i decays at the end of the two 
cascades. Therefore, the missing Et can be expected to 
follow one of the hemisphere axis. It was found to follow 
preferentially the selected bb pair, as can be seen from 
Fig. [14] showing the A(f> distribution between the bb sys- 
tem and the missing Et- Figure [TBI shows the invariant 
mass distribution with Z\</>(Higgs, E^ iss ) < 120°. 

Figures [16] and [17] show the invariant mass distribu- 
tions for the signal and background in the 24-dimensional 
representation with iiia = 190 and 210 GeV/c 2 and with 

E T jet > 100 GeV. 
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Fig. 6: Distribution of Et of the most energetic jet (a), the sub- leading jet (b), and the third most energetic jet (c) 
for the signal with hia = 190 GeV/c 2 and tan/3 = 10 (solid line), for the SUSY background (dashed line) and for the 
tt background (dotted line) in the representation 24. 
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Fig. 7: Distribution of px for the 
b quarks from H/A — > bb (filled 
histogram) and from h — > bb 
(solid line) for the signal with 
niA = 190 GeV/c 2 in the repre- 
sentation 24. 



Fig. 8: Distribution of miss- 
ing Et for the signal with 
m A = 190 GeV/c 2 and tan/3 = 
10 (solid line) in the represen- 
tation 24, for the SUSY back- 
ground (dashed line) and for the 
tt background (dotted line) . 



Fig. 9: Distribution of AR be- 
tween the selected jets for the 
signal in the singlet representa- 
tion with mA = 190 GeV/c 2 . and 
for the tt background. The Z\R 
between the b quarks from the 
Higgs boson decay also shown. 



The selection cuts are summarized below: 

- At least 4 reconstructed jets with Et > 50 GeV 

- One jet with E T > 200 GeV 

- Two jets with E T > 150 GeV 

- Three jets with E T > 100 GeV 

- At least 2 b jets with E T > 50 GeV for the light Higgs 
boson h and at least 2 b jets with Et > 100 GeV for 
the heavy Higgs bosons 

- Missing E T over 200 GeV 

The number of events after each cut and relative cut 
efficiencies are shown in Tables 2] and [5j The QCD back- 
ground was estimated without trigger simulation and with 



full simulation for the efficiency of the E!p iss > 200 GeV 
cut. 



7 Estimation of statistical significance 

The event selection procedure described above yields a vis- 
ible signal in the bb mass as can be seen from Figs.[T2l-[T71 
Therefore the background shape can be determined with 
fits to the side-bands of the peak. Approximately 50% of 
the background level in the signal area is due to the combi- 
natorial background. For the light and heavy Higgs bosons 
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Fig. 10: Invariant mass of the two b jets (filled his- Fig. 11: As in Fig. EDI for 24 representation 

togram) in the singlet representation. The mass with two different E^ et cuts. The mass distribu- 

distribution with correct b jet assignments is tion with correct b jet assignments is shown for 

shown with dashed line. E T jet > 100 GeV. 



Table 4: Cross section times branching fraction, number of events for 10 fb _1 and percentage cut efficiencies (in 
brackets) for the H/A signal with ula = 190 GeV/c 2 and for the SUSY background in the 24 dimensional representation 
with tan/3 — 10. The results for the tt and QCD backgrounds are also shown. For the QCD background the efficiency 
for E^? lss > 200 GeV was calculated with full simulation [48], without trigger simulation. 



rep 24 


Signal 


SUSY 


tt 


QCD 


a x BR (pb) 


5.07 


18.05 


840 


6.29 * 10 4 


Ll+HLT 


36118 (71.2) 


1.51 * 10 s (83.7) 


3.16* 10 s (3.8) 




At least 4 jets, E T > 50 GeV 


23232 (64.3) 


1.28 * 10 b (84.6) 


2.54 * 10 b (80.4) 


1.96 * 10 M 


Et 1SS > 200 GeV 


22871 (98.5) 


96382 (75.4) 


48831 (19.2) 


233949 (0.012) 


E^' cuts 


17315 (75.7) 


53558 (55.6) 


17352 (35.5) 


15334 (6.6) 


At least 2 b jets, E T > 100 GeV 


4663 (26.9) 


4335 (8.1 ) 


3528 (20.3) 


184 (1.2) 


Jet pairing selection 


948 (20.3) 


503 (11.6) 


292 (8.3) 


4.4 (2.4) 



Table 5: Cross section times branching fraction, number of events for 10 fb _1 and prosentual cut efficiencies (in 
brackets) for the h signal with iha = 190 GeV/c 2 and for the SUSY background in the singlet representation with 
tan/3 = 10. The results for the tt are also shown. 



rep I 


Signal 


SUSY 


tt 


a x BR (pb) 


4.13 


14.2 


840 


Ll+HLT 


31782 (64.8) 


119582 (50.4) 


3.16* 10 s (3.8) 


At least 4 jets, Et > 50 GeV 


30682 (96.9) 


96548 (80.7) 


2.54 * 10 5 (80.4) 


e m,ss > 200 Ge y 


20627 (67.2) 


73157 (75.8) 


48831 (19.2) 


E\ et cuts 


15315 (74.2) 


48260 (66.0) 


17352 (35.5) 


At least 2 b jets, E T > 50 GeV 


5494 (35.9) 


5570 (11.5) 


5098 (29.4) 


Jet pairing selection 


1813 (33.0) 


839 (15.1) 


606 (11.9) 


^(Higgs, E£ iss ) < 120° 


866 (69.5) 


254 (51.6) 


235 (61.8) 
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Fig. 12: Distribution of the invariant mass recon- 
structed from the selected b jets for the signal, for 
the total background and for the tt background in 
the singlet representation with iiia = 190 GeV/c 2 . 
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Fig. 13: Distribution of the invariant mass recon- 
structed from the selected b jets for the signal, 
for the total background with Z\R(jj) < 1.5 in the 
singlet representation with iiia = 190 GeV/c 2 . 
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Fig. 14: Distribution of A<f> between the Higgs 
boson reconstructed from the selected b jets and 
missing Ex for reconstructed missing Ex (solid 
histogram) and for missing Ex reconstructed at 
the generator level (solid line) in the singlet rep- 
resentation with iba = 190 GeV/c 2 . 



Fig. 15: Distribution of the invariant mass recon- 
structed from the selected b jets for the signal, 
for the total background with Z\R(jj) < 1.5 and 
zi0(Higgs, E™ 1SS < 120°) in the singlet represen- 
tation with m A = 190 GeV/c 2 . 



with niA = 190 GeV/c 2 in the 24- and 1-dimensional rep- 
resentations the background areas were taken as those 
outside the intervals of 120 < m b g < 250 GeV/c 2 and 
80 < m b 5 < 150 GeV/c 2 , respectively. An interval of 150 < 
m bb < 250 GeV/c 2 was used for m A = 210 GeV/c 2 in 
the 24-dimensional representation. The background dis- 
tributions were obtained fitting a 5 th -order polynomial to 
the sidebands of the mass peaks, the areas outside the 
mass windows given above. The full mass distributions 
were then fitted combining the obtained polynomial and 
a Gaussian distribution. The expected number of back- 
ground events was extracted from extrapolations of the 
background fits to the signal area. The signal events are 
then given as the difference between the measured number 
of events and the expected number of background events 



from the fit. Table [6] shows the mass intervals, the fit- 
ted Higgs boson mass values and the resulting statistical 
significance. Figures [18] - [20] show the fitted mass distri- 
butions for the signal and backgrounds. 

To estimate the statistical significance a likelihood- 
ratio method was used starting from the shapes of the fit- 
ted distributions for the signal+background and for back- 
ground alone. A new "toy experiment" was formed fluctu- 
ating each bin in the distribution according to a Poisson 
distribution with the smeared bin value as a mean. The 
new distribution was refitted with the signal (Gaussian) 
and background (5th order polynomial) shapes. This pro- 
cedure was repeated many (1000) times. For each exper- 
iment, a statistical significance was calculated according 
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Fig. 16: Distribution of the invariant mass recon- 
structed from the selected b jets for the signal, for 
the total background and for the tt background 
in the representation 24 for iiia = 190 GeV/c 2 
with E T > 100 GeV for the tagged b jets. 




rrijj (GeV/c 2 ) 

Fig. 17: Distribution of the invariant mass recon- 
structed from the selected b jets for the signal, for 
the total background and for the tt background 
in the representation 24 for mA = 210 GeV/c 2 
with E T > 100 GeV for the tagged b jets. 




Fig. 18: Fit of a Gaussian plus a 5th degree poly- 
nomial to the bb invariant mass for the signal and 
the total background in the singlet representation 
with m A = 190 GeV/c 2 . 



Fig. 19: Fit of a Gaussian plus a 5th degree poly- 
nomial to the bb invariant mass for the signal and 
the total background in the 24-representation 
with m A = 190 GeV/c 2 . 



Table 6: Higgs boson mass window used for the fits, Higgs boson mass and statistical significance obtained from fits 
with likelihood-ratio method starting from the shapes of the fitted distributions. 





rep24, m A = 190 GeV/c 2 


rep24, m A = 210 GeV/c 2 


repl,m A = 190 GeV/c 2 


Am (GeV/c 2 ) 


150 - 250 


150 - 250 


70 - 150 


niHiggs (GeV/c 2 ) 


193.6±2.9 


224.7±3.7 


118.5±1.5 


Significance 


9.5 


5.6 


21.2 
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Fig. 20: Fit of a Gaussian plus a 5th degree poly- 
nomial to the bb invariant mass for the signal and 
the total background in the 24-representation 
with m A = 210 GeV/c 2 . 



Fig. 21: Distribution of the likelihood-ratio esti- 
mator <7l for q, g —>■ xiH/A — > Xibb in the rep- 
resentation 24 for mA = 210 GeV/c 2 . 



to the following likelihood-ratio estimator: 

o L = ^2ln(L s+B /L B ) (8) 

where -Z/s+b is the maximum likelihood value obtained 
in the signal+background binned likelihood fit, and L-q is 
the maximum likelihood value from the binned background- 
only hypothesis fit. Figure [2X1 shows the distribution of ctl 
for mA = 190 GeV/c 2 in the 24-dimensional representa- 
tion with 10 fb _1 . The distribution has a Gaussian shape 
with width close to unity. The statistical significance can 
be therefore obtained as a mean value of a Gaussian fit 
and is 8.6 for hia = 190 GeV/c 2 . The same procedure 
was performed for iiia = 210 GeV/c 2 . The fitting pro- 
cedure was found to be quite robust and not sensitive 
to the chosen signal region. No systematic uncertainties 
are included. The background shape and the form of the 
function to fit it was assumed to be determined from MC 
studies. The SM backgrounds will be ultimately measured 
from the data. 

8 Methods to measure the backgrounds from 
data 

Methods to measure the SM backgrounds from data have 
been developed in CMS [50] . Missing E T in the QCD 
multi-jet events is mainly due to detector resolution and 
mis-measured jets. Contribution from badly measured jets 
can be suppressed with topological cuts, excluding events 
where the missing Et is parallel to the leading jet. For 
the QCD background, cleaned as much as possible from 
the mis-measured jets, the shape and normalization could 
be obtained measuring the multi-jet events with low Et 
thresholds, where the SUSY contribution is insignificant, 
and extrapolating the Fi^ lss distribution to the area of 
the SUSY searches. Due to large cross section of the QCD 



multi-jet background the contamination from tt events is 
small and could be suppressed demanding that no b-jets 
are present in the sample. 

The large F^ lss values in the tt events with hadronic 
selection are due to detector resolution, mis-measured jets, 
un-identified leptons from W decays, semileptonic decays 
of b and c quarks and hadronic r decays from W — * rv. 
As for the QCD multi-jet background the contribution 
from mis-measured jets can be suppressed with topological 
cuts. The hadronic r decays and the leptons in jets can 
be efficiently identified with the newly developed Particle 
Flow methods [57]. Requirement of two well identified b 
jets and a well measured hadronic top mass may not be, 
however, enough to suppress a contamination from QCD 
multi-jet background for the selection of tt events with 
low Et thresholds. Therefore, the hadronic tt background 
may have to be estimated with Monte-Carlo methods. 



9 Discovery potential 

Figure [22] shows the 5<r-discovery potential for Higgs bo- 
son production in the SUSY cascades g — * qq, q — > Xili 
X2 — * Xih/H/A, h/H/A — > bb in the non-universal 24 
dimensional representation for integrated luminosity of 

10 fhr 1 . Two mass points, m A = 190 and 210 GeV/c 2 
were investigated. The discovery potential was scanned 
by renormalizing the signal and SUSY background cross 
sections as a function of tan/3. The cross sections were 
calculated as described in Section 3. Constant signal and 
background selection efficiency as a function of tan/3 was 
assumed. 

The discovery region extends to masses of ~ 210 GeV/c 2 
around tan/3 = 10 and to ~ 190 GeV/c 2 around tan/3 = 30. 
For ~ 190 GeV/c 2 the low tan/3 region was covered with 
better than 5cr significance down to tan/3 = 2 which was 
the lowest value of tan/3 studied. Since the cross section for 
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H/A production decreases rapidly at Higgs masses around 
~ 280 GeV/c 2 , the discovery region is not expected to 
cover the high Higgs mass values of (mA, tan/3) parameter 
space even with high luminosities. However, large fraction 
of the WMAP preferred relic density region can be ex- 
plored with this channel at the LHC already with 10 fb _1 . 

The systematic uncertanties are not included in the re- 
sult of Fig. [521 The main source of systematic uncertainties 
is due the uncertainty uncertainty of the jet energy scale, 
~ 10% [50], which can affect the efficiency of selecting the 
4 jets and the efficiency of the E^? lss cut. Further sources 
of systematic uncertainty are the b-tagging (~ 5%) and 
the luminosity measurements (~ 5%) [50] . The theoretical 
uncertainty on the tt cross section due to variation of the 
scale and PDF is about 6% [50]. The MC statistics of the 
tt background used in this work corresponds to 2.4 fb _1 , 
which causes an addional uncertainty of about 4% on the 
level of this background after selection. Including the sys- 
tematic uncertainties of similar level has been shown in 
Ref. [50 to degrate the discovery limits in the SUSY pa- 
rameter scape only slighty for integrated luminosities up 
to 10 fb _1 . In particular, the reach for the jet+MET chan- 
nel has been shown to vary by < 5% as a function of mo 
and mi/2 for 10 fb _1 . Similar variation (~ 10 GeV/c 2 ) 
can be expected for the reach of Fig. [22] in iia. 



CMS 

10 fb" 1 




Representation 24 



pp-^H/A+x.,-* bb + X 

WMAP preferred relic density region 



No Systematic Uncertainties 



100 150 200 250 300 350 400 450 500 

m A (GeV/c 2 ) 

Fig. 22: CMS discovery potential at 5a level for q, g — > 
XiH/A — » Xibb in the representation 24 with an inte- 
grated luminosity of 10 fb -1 . No systematic uncertainties 
are included. 



10 Conclusions 

The dark matter allowed regions in the SU(5) GUT rep- 
resentations were studied. From these representations all 
but the singlet lead to non-universal gaugino masses. It 



was demonstrated that combinations of different represen- 
tations can give observed dark matter for otherwise exper- 
imentally allowed parameter values. In the representation 
24 the Higgs signal from the neutralino decay is interest- 
ing and different from the usual universal singlet model. 
Therefore this work especially concentrated on this rep- 
resentation. The preferred dark matter regions for several 
parameter selections were showed and analyzed, including 
the points studied with realistic detector simulations. 

To study the discovery potential at the LHC, the Higgs 
boson, produced in one of the SUSY cascades, was as- 
sumed to decay to b quarks. Event selection was based on 
the requirement of large missing Et, four energetic jets of 
which two tagged b jets. The combinatorial background 
from false jet assignments was reduced by separating the 
jets into two hemispheres. Contributions from the SUSY 
and SM (tt) backgrounds were found to be comparable. 
Visible signals were obtained as a function of the bb invari- 
ant mass. A likelihood-ratio method was used to extract 
the statistical significance. The simulation work showed 
that a 5c-discovery can be expected up to the masses of 
m A = 190 - 210 GeV/c 2 for tan/3 = 5-35. The WMAP pre- 
ferred relic density region, excluding those for tan/3 ~ 8- 
10, can be explored at the LHC in the representation 24 
with the g -> qq, q -> \ 2 q, x~2 -> X\h/H/A, h/H/A -> bb 
cascades in SUSY events with the integrated luminosity 
oflOftr 1 . 
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